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In this review we focus on the idea of establishing connections between the mechanical properties of DNA- 
ligand complexes and the physical chemistry of DNA-ligand interactions. This type of connection is interesting 
because it opens the possibility of performing a robust characterization of such interactions by using only one 
experimental technique: single molecule stretching. Furthermore, it also opens new possibilities in compar¬ 
ing results obtained by very different approaches, in special when comparing single molecule techniques to 
ensemble-averaging techniques. We start the manuscript reviewing important concepts of the DNA mechan¬ 
ics, from the basic mechanical properties to the Worm-Like Chain model. Next we review the basic concepts 
of the physical chemistry of DNA-ligand interactions, revisiting the most important models used to analyze 
the binding data and discussing their binding isotherms. Then, we discuss the basic features of the single 
molecule techniques most used to stretch the DNA-ligand complexes and to obtain force x extension data, 
from which the mechanical properties of the complexes can be determined. We also discuss the character¬ 
istics of the main types of interactions that can occur between DNA and ligands, from covalent binding to 
simple electrostatic driven interactions. Finally, we present a historical survey on the attempts to connect 
mechanics to physical chemistry for DNA-ligand systems, emphasizing a recently developed fitting approach 
useful to connect the persistence length of the DNA-ligand complexes to the physicochemical properties of 
the interaction. Such approach in principle can be used for any type of ligand, from drugs to proteins, even 
if multiple binding modes are present. 
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I. INTRODUCTION 

The DNA molecule is the biological polymer related 
to some of the most important vital processes, from the 
storage and transmission of genetic information to the 
translation of proteins. Its primary structure is usually 
described as two parallel strands with a peculiar chemi¬ 
cal structure based in complementary base-pairs, allow¬ 
ing the replication of the molecule in an unmistakable 
wayi^. The two DNA strands are arranged forming a 
double-helix structure that sets important properties to 
the molecule such as a well-defined negative charge den¬ 
sity and a bending stiffness which places DNA in the class 
of semi-flexible polymers^—. 

Since it stores the genetic information of an organ¬ 
ism, the DNA molecule may be very long in some cases. 
In fact, the human genome has approximately 3 billion 
base pairs, corresponding to a linear contour length of 
the order of 1 meter. If a DNA molecule with this length 
is placed disperse in a water-based solution, its radius 
of gyration will be of the order of 100 /rm£. How can 
a molecule with this size be stored in the nucleus of a 
cell, which has typical dimensions on the order of a few 
micrometers^? The answer lies, at least partially, in the 
mechanical properties of the DNA molecule, which must 
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be unique to allow such a condensation. In vivo , this 
process usually occurs mediated by the interaction of the 
DNA molecule with ligands, especially (but not exclu¬ 
sively) histone proteins. Furthermore, from molecular 
biology it is known that other important intracellular 
processes such as cell division and protein binding also 
depend on the DNA topology, which in turn, depends 
on the mechanical properties of the DNA molecule^. 
DNA topology can be strategically changed during these 
processes by the action of enzymes such as the topoiso- 
merases, allowing their occurrence efficiently^!. 

Like the proteins and enzymes exemplified above, 
many drugs are capable to interact with DNA, modify¬ 
ing its mechanical properties with biological implications 
in vivo. Cancer chemotherapy, for instance, is a field 
in which the details about DNA interactions with drugs 
are important. In fact, some classes of drugs such as 
the anthracyclines and the platinum-based compounds 
exhibit a strong affinity to interact with the DNA of can¬ 
cer cells. When these drugs bind to DNA they can in¬ 
hibit the replication process, thus stopping the tumor 
growthiiii- . On the other hand, gene therapy is another 
field of medical sciences in which this kind of knowledge 
is also importan t 13 ’ 14 . In these therapies, DNA molecules 
are usually transported from outside to inside living cells 
in order to replace defective genes, thus correcting cell 
malfunctions. One approach to accomplish this transport 
easily, for example, is condensing the DNA molecule by 
using cationic ligand s 15 ' 16 . 
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In summary, all the examples discussed above show 
the importance in studying and understanding the details 
behind DNA interactions with ligands. In fact, many 
researchers of varied areas such as physics, chemistry, 
biology, medicine, pharmacy, engineering, etc have paid 
attention to this topic along the past 20 years, with a fast 
increase of the number of publications and citations^. 

In this manuscript we review important topics of the 
field “DNA-ligand interactions”, emphasizing in how one 
can connect the changes of the mechanical properties of 
the DNA induced by the binding ligand to the physic¬ 
ochemical information of such interaction. This type of 
connection is interesting because it allows one to perform 
a robust characterization of the interaction both from the 
mechanical and physicochemical point of view by using 
only one experimental technique: single molecule stretch¬ 
ing experiments. To discuss such connection, firstly in 
Section 2 we discuss the basic DNA mechanics, revisit¬ 
ing the main concepts and approaches used in the field. 
In particular, we revisit the Worm-Like Chain (WLC) 
model, the standard one used to describe bare DNA me¬ 
chanics and to investigate the changes of the DNA me¬ 
chanical properties when interacting with a binding lig¬ 
and. Then, in Section 3 we discuss the physical chem¬ 
istry of DNA-ligand interactions, emphasizing the chem¬ 
ical equilibrium states which can usually be described by 
a binding isotherm. We revisit the most important mod¬ 
els used in the field, discussing their strong points and 
limitations. In Section 4 we discuss briefly the experi¬ 
mental techniques most used to perform single molecule 
experiments, emphasizing the key features of each one. 
In Section 5 we present and discuss the main types of in¬ 
teractions that occur between DNA and ligands: interca¬ 
lation, covalent binding, electrostatic driven interactions 
and groove binding. Finally, in Section 6 we present and 
discuss the main topic of this review: the approaches on 
how one can connect mechanics to physical chemistry. 
The final conclusions are presented in Section 7. 


II. DNA MECHANICS 

During the last decades, DNA mechanics has become 
a very well studied topic especially due to the advent of 
single molecule techniques. Such techniques allow one to 
manipulate and stretch individual DNA molecules, giving 
access to mechanical information contained in the “force 
x extension” curves. Before single molecule techniques, 
such type of information was somewhat difficult to be 
accessed by ensemble-averaging techniques. 

From the mid-90s some theoretical models were formu¬ 
lated in order to explain the mechanical behavior of DNA 
molecules. In particular, most of these models attempt to 
give a theoretical expression for the “force x extension” 
curve based on key mechanical parameters such as the 
linear contour length of the polymer chain and the DNA 
bending stiffness, which can be conveniently represented 
by its persistence length. 


The contour length is the most basic mechanical prop¬ 
erty of a polymer chain: it is simply the length of the 
chain measured along its contour, which is proportional 
to the number of monomers. The persistence length, oth¬ 
erwise, is the correlation length of the polymer chain and 
thus gives information about the bending stiffness of the 
polymer. In the case of DNA molecule, the persistence 
length has basically two components: the intrinsic and 
the electrostatic one. The first component is related to 
the bending rigidity due to the molecule composition it¬ 
self, while the second one is due to the negative charge 
distribution along the double-helix=2r 2 i. Since these two 
components are usually present in most relevant situa¬ 
tions, the models in general represent the persistence 
length by its effective value, which takes into account 
the two contributions. Following most authors, in this 
manuscript we will call the effective persistence length 
only by persistence length. In water-based solutions un¬ 
der nearly physiological conditions (pH = 7.4, [NaCl] = 
150 mM), a disperse bare DNA molecule is classified as a 
semi-flexible (or semi-rigid) polymer due to its interme¬ 
diate value of the bending stiffness, which corresponds to 
a persistence length A ~ 50 nn& 22 r— . 

In the following section we present the most relevant 
model used for studying the mechanics of DNA molecule: 
the Worm-Like Chain (WLC) model. In this manuscript 
we do not intend to review other mechanical models or 
present an historical survey on this specific subject, since 
today the WLC model is recognized as the standard one 
to study DNA mechanics. 

A. Worm-Like Chain Model (WLC) 

The Worm-Like Chain (WLC) is a model derived from 
polymer physics, and has become in the past years the 
standard one in analyzing DNA stretching experiments. 
To introduce this model, let us firstly assume that the 
polymer itself is a chain formed by rigid rods with lengths 
b, connected by freely-rotating vertices. Let us call 0i the 
angle between the rods i and i + 1. The WLC model 
is then defined by assigning an harmonic bending energy 
function to the angle formed between the two rods^— , 

m) = Y b °i’ (!) 

taking the continuum limit with b —t 0. The constant k 
is the effective elastic bending stiffness of the chain. 

In the continuum limit, Eq. |T]can be used to write the 
total bending energy of the chaiiA2i, 

E = \ (2) 

where C is the local curvature at each point, ds is a length 
element along the polymer, and L is the contour length 
of the polymer chain. 
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The parameter k is directly related to the polymer per¬ 
sistence length A by 


exact numerical solution of the WLC models, which was 
done perturbatively. The resulting expression reads 


A = 


K 


( 3 ) 


where k B is Boltzmann’s constant and T is the absolute 
temperature. 

Equation [2] can be used to deduce the behavior of 
the force as a function of the polymer extension as one 
stretches it. This analysis can be performed numerically 
or analytically using appropriate approximations^. In 
1995, Marko and Siggia solved the model analytically, 
obtaining an approximate expression for the force as a 
function of the polymer extension which has become the 
most used to analyze DNA stretching experiments in the 
entropic low-force regime (F < 5 pN)^2. Their result is 


F = 


k B T 

A 


z 

L 


1 

4 (!-f ) 2 


1 

4 


( 4 ) 


where F is the force and z is the end-to-end distance 
(extension) of the DNA molecule. 

Despite its renowned utility, this expression is still an 
approximation, diverging at z = L. Moreover, Eq. [|] 
describes well only the entropic regime of the polymer, 
which is valid for stretching forces typically below ~ 5 
pN. In this regime the applied forces are sufficiently small 
such that they can change only the polymer conformation 
in solution, i. e., its entropy. 

Also in 1995, Odjik proposed a different approach that 
accounts for higher forces, in which enthalpic effects 
start to became relevant for the polymer mechanics^. 
The enthalpic regime is defined as the regime in which 
the stretching forces became large enough to distort the 
DNA primary structure and eventually to break chemi¬ 
cal bonds. Such effect can be accounted by introducing 
an enthalpic mechanical parameter to describe the poly¬ 
mer deformation: the stretch modulus S. The analytical 
expression proposed by Odjik reads^ 


z = L 

Observe that the stretch modulus S has units of force. 
Taking the limit S —> oo and inverting the above equa¬ 
tion (isolating F), we found an equation similar to the 
Marko-Siggia expression (Eq. [4| if z ~ L, i. e., neglect¬ 
ing very small forces. Thus, observe that a polymer in 
the entropic regime can be interpreted as a polymer that 
has a stretch modulus S very large, i. e., that resists 
deformations on its chemical structure. 

In 1999, Bouchiat et al. proposed another solution of 
the WLC model in the entropic regime. Their approach 
consists in adding six terms to Eq. [I] in order to improve 
its accuracy^. These terms were determined by com¬ 
paring the results predicted by Eq. [J]to results from an 


1 l k ^ T . F 

2 \l AF + S 


( 5 ) 


F = 


k B T 

A 


z 

L 



where the a^s are constants numerically determined. 

In addition to the models discussed above, important 
contributions to the elucidation of many peculiarities of 
DNA mechanics were given by the groups of A. Volo- 
godskii, M. D. Frank-Kamenetskii, H. E. Gaub, M. C. 
Williams, V. Croquette, F. Ritort, C. Bustamante and 
others, especially concerning the bending of small DNA 
fragments, strong bending and fluctuations in the double¬ 
helix, dependence of DNA rigidity on the temperature 
and base sequence, DNA twist, overstretching transition, 
DNA hairpins, etci 21 i 22 ' 28 i 30 ~ — . 


III. PHYSICAL CHEMISTRY OF DNA-LIGAND 
INTERACTIONS 

The study of the physical chemistry of DNA-ligand in¬ 
teractions consists in two different sub-fields: the chemi¬ 
cal equilibrium of the interaction and the kinetics of the 
interaction. Consider the system of interest (DNA -f- lig¬ 
and molecules in solution) as composed by two different 
partitions where the ligand molecules can stay: the DNA 
(bound ligand molecules) and the solution (free ligand 
molecules). The chemical equilibrium is achieved when 
the average number of molecules in the partitions remains 
constant in time. The kinetics of the interaction, other¬ 
wise, describes the changes that occur between the initial 
incubation and the final equilibrium state. 

In this manuscript we emphasize the physical chem¬ 
istry of the chemical equilibrium, since the equilibrium 
states can be represented by a binding isotherm that can 
be linked to the changes of the mechanical properties 
of DNA-ligand complexes. Below we discuss the most 
relevant models that attempt to describe the chemical 
equilibrium of DNA-ligand interactions. Some studies on 
the kinetics of such interactions were performed by the 
groups of M. C. Williams, D. Anselmetti, D. M. Crothers 
and others^ - —. 


A. The general problem 

Consider two molecules A and B associating in solu¬ 
tion to result in a molecule C. This mechanism can be 
represented by the chemical reaction 


A + B U C, (7) 

K d 

where Ki and Kd are, respectively, the equilibrium in¬ 
trinsic binding constants of association and dissociation. 
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They are also known as thermodynamic constants or 
macroscopic constants. Observe that AT* represents the 
association reaction, where the reagents A and B asso¬ 
ciate to result in the compound C, while Kd represents 
the dissociation reaction, i. e., the reverse reaction in 
which C dissociate in the original reagents A and B. 

These constants are defined in term of the molar con¬ 
centrations of the involved substances, 


Ki = 


[C] 

[A][By 


( 8 ) 


only for very small ligand molecules which occupy only 
one DNA base-pair when bound, which is not the case 
for most ligand molecules, (b) It supposes that previous 
bound ligand molecules do not interfere in the binding 
mechanism of the subsequent ones, i. e., the interaction 
is non-cooperative. 

The first simplification can be bypassed by introducing 
the parameter r moj , the bound ligand fraction at satura¬ 
tion, i. e. the maximum value of the bound ligand frac¬ 
tion r. Observe that the inverse of r max is the mean num¬ 
ber of base pairs occupied by each bound ligand molecule 
N = 1 j t max ■ The corrected binding isotherm then reads 


and 


K d 


[MB] = 
[C] 


(9) 


Note that in these last two equations, [A] is the molar 
concentration (1 M = 1 mol/liter) of the compound X. 
Also observe that AT,: has units of M -1 , while Kd has 
units of M. 


B. Scatchard model 

This is the simplest model that describes the chemical 
equilibrium of the DNA molecule with ligands in solu¬ 
tion. Let us firstly adapt the previous notation for the 
specific case of DNA-ligand interactions. Call [A] = Cf 
the concentration of free ligands solution and [ C] = C b 
the concentration of ligands bound to DNA (result of 
the reaction). Suppose firstly that each ligand molecule 
occupies only one base pair of the DNA when bound. 
Consequently, the concentration of free linkable sites in 
the DNA molecule can be written as [ B] = Cb P - Cb, 
where Cb P is the concentration of DNA base pairs, which 
is a constant. 

Substituting these definitions in Eq. [HI one has 


Ki = 


C b 

Cf(C bp -C b )- 


( 10 ) 


Now we introduce the bound ligand fraction r, 


r = 


Cb_ 

c bp ' 


such that Eq. [TO] can be rewritten as 


( 11 ) 


KjCf 
1 + KiCf ’ 


( 12 ) 


which is known as the Scatchard binding isotherm, pro¬ 
posed originally in 1949^. 

Despite its didactic utility, the Scatchard binding 
isotherm has two important simplifications: (a) It is valid 


I rnax K j C f 
1 + K,C f ’ 


(13) 


C. Hill model 

The Hill binding isotherm was originally proposed by 
A. V. Hill in 1910 to describe the binding of oxygen to 
hemoglobin inside red blood cellsSS,. 

Basically the model introduces the Hill exponent n, a 
cooperativity parameter which is a lower bound for the 
number of cooperating ligand molecules involved in the 
reactio n 63 : 64 . The binding isotherm reads 

_ r max {KiCf) n 

1 + {KiCf) n ' 1 ; 

The apparent binding association constant of the re¬ 
action is defined as Ka = AT". Observe that if n > 1, 
the interaction is positively cooperative, i. e., a bound 
ligand molecule increases the apparent affinity of DNA 
for subsequent ligand binding. If n < 1, otherwise, the 
interaction is negatively cooperative and a bound ligand 
molecule decreases the apparent affinity of DNA for sub¬ 
sequent ligand binding. If n = 1, the interaction is non- 
cooperative and the affinity is independent of the number 
of previously bound ligand molecules. 

The Hill binding isotherm has achieved a particu¬ 
lar success to describe positively cooperative “none-or- 
all” processes in > 1), in which the cooperating ligand 
molecules bound practically simultaneously to the bound 
site forming a bound cluste r 63 ! 65 . On the other hand, 
when n = 1 the Hill isotherm reduces to the Scatchard 
one and is therefore able to describe individual binding 
of ligand molecules^. Finally, to the best of our knowl¬ 
edge there is no report in the literature of a negatively 
cooperative DNA-ligand interaction described by a Hill 
binding isotherm. 


D. Neighbor exclusion model (NEM) 

This model was proposed in 1974 by McGhee and von 
Hippel with the purpose of analyzing in detail the neigh¬ 
bor exclusion effects due to large ligand molecules that 
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occupy more than one DNA base-pair— >— The authors 
have accounted for the ligand size by introducing the ex¬ 
clusion parameter N, the number of base-pairs that a lig¬ 
and molecule effectively occupies when binding to DNA. 
This parameter was cited earlier in connection to the sat¬ 
urated bound ligand fraction, N = l/r max . 

The model has a non-cooperative and a cooperative 
version, but the last one has been used only in a few 
work s 68 ' 69 to analyze experimental data because the 
binding isotherm is somewhat intricate. 

The non-cooperative binding isotherm reads 


Cf 


= Ki( 1 - Nr) 


1 — Nr 
1 - {N - 1)7 


-i JV-l 


and the cooperative binding isotherm reads 


(15) 


Cf 


with 


= Ki( 1 - Nr) 


(2w —1)(1 — Nr) + r-R 
2(w - 1)(1 - Nr) 

1 - (N + l)r + R 
2(1 - Nr) 


N-l 


U6) 


R = a s /[1 — (N + 1 )r] 2 + 4wr(l — Nr). (17) 

Here ui is the cooperativity parameter. For u> smaller, 
equal, or larger than unity, one has negative, non- 
cooperative, or positive cooperativity, respectively. 

The major advantage of this model is to treat in more 
detail the effects related to the ligand size. This feature 
is particularly important in the analysis of DNA inter¬ 
actions with intercalators, a class of ligands in which 
neighbor-exclusion effects is extremely important— - —. 
In fact, NEM has become in the past years the stan¬ 
dard binding isotherm used to analyze DNA interactions 
with intercalatorsSU—. 


IV. SINGLE MOLECULE EXPERIMENTAL METHODS 

In this section we briefly discuss the single molecule 
experimental techniques commonly used to measure the 
mechanical properties of the DNA-ligand complexes. 

The main advantage of single molecule techniques is 
the possibility to study a particular DNA molecule free 
from the influence of other molecules in the sample. Sin¬ 
gle molecule stretching experiments such as those per¬ 
formed with optical or magnetic tweezers usually give in¬ 
sights on the global (long length scale) mechanical prop¬ 
erties of individual DNA molecules. In fact, mechanical 
parameters such as the persistence and contour lengths 
and the stretch modulus can be extracted by analyzing 
the force x extension curves of the complexes, which can 
be obtained in single molecule approaches. 

Useful reviews which discuss and compare single 
molecule techniques can be found in the literature^ - — . 


A. Optical tweezers 

Since the seminal works of Ashkin and 
collaborators^^, optical trapping and manipula¬ 
tion have found various applications in many areas of 
science such as physics, biology and chemistry. Today, 
the most common optical tweezers are mounted by 
focusing a laser beam with a microscope objective of 
large numerical aperture. This apparatus can trap small 
dielectric objects near the lens focus, being a powerful 
tool to manipulate beads, particles and biological 
systems with typical sizes in the micrometer rang o 82 ' 83 . 
The typical forces obtained with this apparatus are 
between 0.1 - 400 piconewtons, which are in the range of 
many biological forces such as the entropic and enthalpic 
forces on biopolymers and molecular motors. For an 
introductory review about the basic theory and features 
of optical tweezers, see ref.—. Other useful reviews on 
instrumentation and recent advances on the technique 
can also be found in the literature^ - —. 

To perform precise quantitative measurements with 
optical tweezers, size-calibrated dielectric beads have be¬ 
come the standard objects to be captured because of their 
perfect symmetry which facilitates trap calibration and 
position detection. A dielectric bead trapped in an opti¬ 
cal tweezers is an overdamped Brownian harmonic oscil¬ 
lator, such that the optical trap can be characterized by 
its trap stiffness n which depends on the bead size and 
refractive inde:»— 

In the last decades, optical tweezers have been largely 
used to study the mechanical properties of DNA/RNA 
molecules and their complexes formed with drugs or 
proteins. Useful reviews on this subject can be found 
in the literature^ -92 . Basically, the classic experiment 
consists in attaching one end of the DNA molecule 
to a polystyrene or silica bead and the other end to 
a substrate (a microscope coverslip or a second bead 
attached to a micropipette, for example). The opti¬ 
cal tweezers is then used to trap the bead and so the 
DNA molecule can be manipulated and stretched by 
moving the laser beam or the microscope stage. The 
force as a function of extension can be measured as one 
stretches the DNA molecule. To perform this task, one 
needs to detect the bead position and to calibrate the 
tweezers (determine the trap stiffness k). There are 
many techniques which can be used to perform this kind 
of measurement, such as dynamic light scatterin g 93 ' 94 , 
back-focal plane interferometry^,, statistics of thermal 
fluctuations^, simple videomicroscop y 64 ' 97 ' 98 , calibra¬ 
tion using hydrodynamic drag forces^ or by using other 
types of detectors— For a recent review on measuring 
with optical tweezers, see ref—. 

Figure Q] shows a typical force x extension curve of a 
single bare A-DNA molecule (~ 48,500 base-pairs) ob¬ 
tained by performing a DNA stretching experiment in 
the entropic regime with optical tweezers. The trap cali¬ 
bration and the bead position detection were performed 
in this case by using videomicroscopy^i, and the solid 
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DNA extension (|rm) 


FIG. 1. Force x extension curve of a bare A-DNA molecule 
in the entropic regime. Circles', experimental data obtained 
with optical tweezers; Solid line: a fitting to the Marko-Siggia 
Worm-Like Chain (WLC) model (Eq. [4|. For this particular 
DNA molecule we have found from the fitting A = (50 ± 2) 
nm and L = (15.6 ± 0.1) /jrn. 

line corresponds to a fitting using the Marko-Siggia WLC 
model (Eq. |4]). From the fitting one can promptly de¬ 
termine the persistence and contour lengths of the DNA 
molecule, obtaining for this particular curve A = (50 ± 
2) nm and L = (15.6 ± 0.1) /im. 

B. Magnetic tweezers 

The idea behind magnetic tweezers is very similar to 
its optical analogue, the main difference is that in this 
case the forces are exerted by an external magnetic field 
applied around the sample. Paramagnetic beads are used 
instead of dielectric ones in order to be manipulated with 
the magnetic held. The typical forces obtained are of 
the order of hundredths of piconewtons to hundreds of 
piconewtons. 

Basically, the force applied on the paramagnetic beads 
can be written as 

f = (18) 

where ~jt is the magnetic dipole moment induced in the 
bead and ~e5 is the applied magnetic held. Observe that 
for moderate magnetic helds one has oc and the 
resulting force is proportional to the gradient of the held 
intensity. 

Reviews on the basic and advanced features of mag¬ 
netic tweezers can be found in the Iiterature22r-i21. 

An advantage of this technique in relation to opti¬ 
cal tweezers is its convenience to apply torques on the 
magnetic beads by rotating the external magnetic held, 
which allows one to rotate DNA molecules and therefore 


to study quantities such as the torsional rigidity and the 
degree of supercoilin g 22 ' — 2 r225 . . These quantities are 

also mechanical properties important to some biological 
processes in which the double-helix must be unwound, 
such as in DNA replication. Another advantage of the 
magnetic tweezers is its convenience to perform constant- 
force experiments, working as a force-clanrp trap (it is 
just a matter of choosing the adequate magnetic held - 
see Eq. HID. Constant-force experiments can be per¬ 
formed with optical tweezers only using non-conventional 
(and more intricate) approaches such as by using a force- 
feedback electronics or working in anharmonic regions of 
the optical potentia l 106 ’ 107 . Among the disadvantages of 
using magnetic tweezers, one can cite the hysteresis of 
the magnetic held and heat generation around the sam¬ 
ple if the held is produced by current distributions, aside 
the more intricate calibration of the apparatus and its 
restriction to applications with magnetic materials. A 
recent work by Neuman and Nagy provides a detailed 
comparison between optical and magnetic tweezers, and 
also atomic force microscopy^. 


C. Atomic force microscopy (AFM) 

Atomic force microscopy (AFM) is another impor¬ 
tant tool in single-molecule studies of DNA-ligand 
interaction s 108 ’ 109 . In the last 20 years, a number of dif¬ 
ferent protocols have been developed in order to deposit 
DNA molecules on a hat surface and to image them reli¬ 
ably and reproducibly. Today, the standard surfaces used 
to deposit the DNA molecules are mica substrates and 
less often silicon substrates, because of their low rugosity. 
A number of buffer solutions containing divalent cations 
(such as Mg 2+ , Ni 2+ , Mn 2+ , Co 2+ and Ca 2+ ) have been 
used to enhance the DNA adsorption onto the substrate, 
which is otherwise poor. Moreover, divalent cations also 
allow the polymer chain to equilibrate on the hat 2D ad¬ 
sorbing surface, preventing chain kinetic trapping, which 
must be avoided in order to study equilibrium proper¬ 
ties of the adsorbed DNA or DNA-ligand complex*!^. 
Once adsorbed, DNA molecules and DNA-ligand com¬ 
plexes can be imaged using the AFM usually operating 
in the tapping mode, which minimizes possible damages 
to the sample due to the tip-surface interactions during 
the scanning. The images obtained are topographical 
maps which associate a certain height to each point on 
the sample. By analyzing these images, several DNA 
statistical parameters such as the mean contour length, 
the persistence length and bending angles can be esti¬ 
mated. As discussed by Rivetti et aL —, this analysis 
can be performed, for example, by measuring the mean- 
squared end-to-end distance < R 2 > of the polymer. In 
fact, statistical mechanics of polymers predicts that for 
2D worm-like chains (which is the case of deposited DNA 
molecules), < R 2 > is given by 















Extracting physical chemistry from mechanics: a new approach... 


7 


< R 2 >= AAL 


1 - 


2 A / l \ 

TV 1 " 6 "”) 


(19) 


By using this equation it is possible determine the per¬ 
sistence length A by measuring < R 2 > and the contour 
length L for the deposited DNA molecules. 

On the other hand, the visualization of DNA conden¬ 
sates formed with polycations with the AFM technique is 
nowadays a routine in many laboratories. The morphol¬ 
ogy of these condensed DNA complexes seen in the AFM 
images are as well as clear and well-defined as the images 
produced by other kinds of microscopy techniques such as 
electron microscopy (EM). The structure of DNA-protein 
complexes has also been a target of a number of stud¬ 
ies. Within the limits of the AFM technique are, for 
instance, the visualization of sharp kinks and cross-links 
introduced in the DNA molecule by histone-like proteins, 
the structure of nucleosome particles, the determination 
of protein binding-sites and more recently the determi¬ 
nation of protein association constants to DNAiii. 

Finally, besides being a powerful tool for visualizing 
single molecules, the AFM apparatus can also be used 
to perform force spectroscopy in solution like optical or 
magnetic tweezers, allowing one to determine the force x 
extension curves of the DNA-ligand complexe s 36 ! 112 ’ 113 . 

Specific reviews about the application of the AFM 
technique to single molecule studies can be found in the 
literature 1 ^’ 1 

As a final remark, with the improvement of fluorescent- 
based optical technology along the last decades, fluores¬ 
cence microscopy has also became another important tool 
to visualize DNA structure, conformation changes and 
interactions with ligands at single molecule leveli^ 123 . 
The technique can be used as complementary to AFM, 
with the advantage that one does not need to deposit the 
molecules in a substrate. 


V. DNA-LIGAND INTERACTIONS 

DNA can interact with ligands in many different ways, 
from covalent binding to simple electrostatic driven inter¬ 
actions. Here we describe briefly the most relevant types 
of interactions, discussing the main features of each one. 


A. Covalent ligands 

The covalent binding of drugs to DNA is usually ir¬ 
reversible and completely inhibit DNA processes. The 
platinum-based compounds are examples of drugs which 
can interact with DNA by covalent binding-L 1 -^. Cis- 
platin and its related compounds carboplatin and ox- 
aliplatin are antitumor platinum-based molecules usu¬ 
ally used in cancer chemotherapy. The action of these 
complexes as anticancer drugs consists in damaging the 
DNA molecule with adducts that form various types of 


crosslinks, which introduce strong structural perturba¬ 
tions and impede DNA replicatio n 126 ^ . ? . The clinical 
use of these complexes, however, is limited due to their 
several side effects and the development of drug resis¬ 
tance. 

Another example of covalent binding is found in the 
interaction of drugs from the class of furocoumarins (pso¬ 
ralen, angelicin, etc) with DNA when one illuminates 
the complex with ultraviolet-A (UVA) lights. Pso¬ 
ralen is a well-known drug used in the treatment of skin 
diseases like psoriasis, vitiligo, and some other kinds of 
dermatitis^. The most common therapy is called PUVA 
(psoralen followed by UVA light), which consists in tak¬ 
ing a medicine containing psoralen and exposing the pa¬ 
tient to UVA light. The drug effectively increases the 
skin sensitivity to UVA and the skin melanin leve l 130 ’ 131 . 
It is well established in the literature that when a DNA- 
psoralen complex is illuminated with UVA light, the drug 
molecules absorb photons and form covalent bonds pref¬ 
erentially with the thymine s 128 ! 132 . When there is no illu¬ 
mination at the sample, however, psoralen interacts with 
DNA by intercalative binding - see next section. The ef¬ 
fects of covalent binding on the mechanical properties of 
the DNA-psoralen complexes were recently studie d 94 ’ 133 . 
In particular, it was shown that the contour and persis¬ 
tence lengths of the complexes depend on the psoralen 
concentration and on the exposure time to UVA lights. 


B. Intercalators 


Intercalative binding is one of the most common in¬ 
teractions between DNA and ligands, and was firstly de¬ 
scribed by L. S. Lerrnan in 196 1 134 ’ 135 . It is characterized 
by the insertion of a flat aromatic molecule between two 
adjacent DNA base pairs. The complex is thought to be 
stabilized by the stacking interactions between the lig¬ 
and and the DNA basest. Intercalators also introduce 
strong structural perturbations on the double-helix struc¬ 
ture. To accommodate the intercalated molecules, there 
is an increase in the DNA contour length, which is ac¬ 
companied by an unwinding of the double-helix by a cer¬ 
tain angle per intercalated molecul o 56 i 70 ~ - 3,137 . Dauno- 
mycin, doxorubicin and ethidium bromide (EtBr) are 
classic examples of drugs which intercalate in the DNA 
molecule and can modify its elasticity depending on the 
drug concentration. Daunomycin and doxorubicin are 
anthracycline antibiotics used in the treatment of vari¬ 
ous cancers such as some types of leukemias, sarcomas, 
lymphomas, myelomas, neuroblastomas, as well as can¬ 
cers in the breast, head, ovary, pancreas, prostate, stom¬ 
ach, liver, lung and others. They inhibit DNA replica¬ 
tion and transcription when intercalating, impeding cell 
duplication 73 . Ethidium bromide (EtBr) is commonly 
used as a fluorescent stain for identifying and visualizing 
nucleic acid bands in electrophoresis and in other meth¬ 
ods of nucleic acid separation. Other known intercalators 
are the DNA fluorescent stains acridine orange, metliy- 
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lene bluo^ and diaminobenzidine^. More examples and 
specific reviews on the basic properties of intercalators 
can be found in the literature^ - —. 

Many aspects of the DNA-daunomycin interaction, 
such as kinetics, self association and equilibrium bind¬ 
ing were studied by J. B. Chaires, D. M. Crothers and 
coworkers in the 80’ o 60 ' 73 ' 137 ' 139 . On the other hand, the 
DNA-EtBr interaction was characterized in various as¬ 
pects by many authors, but even today one can found 
somewhat contradictory results about the mechanical be¬ 
havior of such complexe s 23 ' 56 ! 75 ' 77 ! 140- -^, and also for 
different complexes formed between DNA and other in¬ 
tercalators, especially when comparing results obtained 
from different experimental technique s 56 ' 75-77 i 115 1 140 ~ i^. 

Recently our group has studied in detail the changes 
in the persistence and contour lengths of DNA complexes 
formed with various intercalating molecules^ -77 ' 94 ! 133 , 
by using optical tweezers in a very low force regime (F < 
2 pN). We reported an abrupt structural transition in 
the persistence length due to drug intercalation, which is 
probably related to a partial denaturation of the DNA 
molecule due to the pulling force used to stretch the 
complexe o 76 i 77 ' 133 i 147 . The contour length, otherwise, 
does not present such a transition, increasing monotoni- 
cally with drug concentration until saturation. 


C. Electrostatic driven interactions 

Since the DNA molecule has a high negative charge 
density in aqueous solution due to its phosphates (2 el¬ 
ementary charges per each 3.4 A along the DNA axis), 
it strongly interacts with itself (i. e., different DNA seg¬ 
ments strongly repel each other thus promoting the chain 
swell) as well as with positively charged ligands such as 
ions and macro-ions, especially multivalent cations!^. 

DNA condensation due to multivalent cations is a clas¬ 
sic example which shows the importance of electrostatic 
driven interactions in DNA solutionsi^hlsl, In this pro¬ 
cess the multivalent cations binds along the DNA double¬ 
helix, and the strong positional correlations between 
them start to play a role and promotes a coil-globule 
transition: the DNA molecule folds onto itself^ ^M with 
a high increase in the local DNA segment density at 
the level of both mono-molecular collapse or in a multi- 
molecular aggregation. 

Some models concerning electrostatic interactions be¬ 
tween DNA and ligands were proposed along the last 
decade s 19 ' 155 ' 156 . In fact, there are different hypotheses 
to explain the DNA bending mechanism by multivalent 
cations, including a purely electrostatic model by Rouz- 
ina and Bloomfield^^ and an asymmetrical phosphate 
neutralization model by Manning^. According to Rouz- 
ina and Bloomfield, a multivalent cation binds to the en¬ 
trance of the DNA major groove, between the two phos¬ 
phate strands, electrostatically repelling sodium counte¬ 
rions from the neighboring phosphates. The unscreened 
phosphates on both strands are strongly attracted to the 


groove-bound cation. This binding leads to groove clo¬ 
sure, accompanied by DNA bending towards the cationic 
ligand^. Differently, Manning proposes that the sta¬ 
ble double-helix structure of DNA represents an equi¬ 
librium between stretching forces (caused by interphos¬ 
phates repulsion) and compressive forces (caused by at¬ 
tractive interaction between nucleotides). This analysis 
suggests that significant local interphosphate stretching 
forces balance compressive forces within DNA and that 
these stretching forces can drive DNA deformation when 
phosphates charge are locally neutralized. 

These two approaches predict a reduction of the per¬ 
sistence length as the concentration of bound cations in¬ 
creases. In fact, the model proposed by Rouzina and 
Bloomfield predicts that the effective persistence length 
Ae of the DNA-ligand complex is given by 


1 _ 1 Nr 
Ae A\ A 2 ’ 


( 20 ) 


where A\ is the bare DNA persistence length (when no 
ligands are bound r = 0) and 1/A\ + 1/A 2 is the in¬ 
verse persistence length of a DNA saturated with ligands 
(which occurs when r = r max = 1/iV). Observe that 
here N is the exclusion parameter of the ligand and r = 
Cb / Cb p is the ratio between the bound ligand concentra¬ 
tion and the DNA base-pair concentration, as introduced 
in Section III. 

The model developed by Manning, on the other hand, 
predicts that the effective persistence length Ae of the 
DNA-ligand complex (charge neutralized DNA) is related 
to the original persistence length Aq (fully charged DNA) 
by the equatio n 19 ' 156 


Ae = 


irR 2 


pA 0 


1 3/2 


2(£ - 1) - ln(re&) 


( 21 ) 


where R is the radius of the double helix, f3 is the Bjerrum 
length (distance between two unit charges in pure solvent 
- no other ions - at which the electrostatic energy is AbT), 
1/re is a measure of the extent of the ion cloud around the 
object, b is the average axial distance between phosphates 
(0.17 mn) and £ = (3/b is a measure of the axial charge 
density of the DNAi— 

This model predicts, for example, that for 30% of neu¬ 
tralized charge, the effective persistence length is Ae = 
33.2 nm. For 60% of neutralized charge, Ae = 11.1 nrn 
and for 100% of neutralized charge, Ae = 7 nm—i. 


D. Major and minor groove ligands 

Most drugs that interact electrostatically with DNA 
usually exhibits a preference to the major or minor groove 
floor of the double-helix. Many minor groove ligands are 
known by their antitumor and antibiotic functions. This 
kind of interaction is usually characterized by a com¬ 
bination of electrostatic, van der Waals and hydrogen 
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bonds. Examples of minor groove ligands are the anti¬ 
cancer compound distamycin A, the antibiotics netropsin 
and berenil, and the fluorescent stain DAPI. These drugs 
usually form reversible complexes with DNA, preferen¬ 
tially binding at AT base pairs sequences. They also in¬ 
duce elasticity changes on the DNA molecule, stabilizing 
the double-helix structure^. An extensive review on the 
DNA minor groove complexes can be found in ref. 1 ' 37 . 

On the other hand, major groove binding is also a 
kind of interaction usually characterized by electrostatic 
binding^. a-Helical (Ac-(Leu-Ala-Arg-Leu)3-NH linker) 
is a peptide which interacts with DNA via major groove 
binding^. Other known examples are the intercalator 
and major groove ligand YO^, the bis-intercalator and 
major groove ligands YOYO and ditercaliniu m 145 1 158 1 159 
and the anticancer drug neocarzinostatini^. More ex¬ 
amples and a discussion on the main characteristics of 
major groove binding ligand can be found in a recent 

1 59 

reviewia*. 


E. Ligands with multiple binding modes 

There are many ligands which can interact to DNA by 
different binding modes, depending on factors such as the 
properties of the surrounding buffer solution, the DNA 
base-pair sequence, external conditions such as sample il¬ 
lumination, etc. Some examples were already cited in the 
last sections. In some cases the ligand has distinct por¬ 
tions which interact to DNA by different modes. In other 
cases there is only a single binding mode for the entire 
ligand molecule, which can be changed upon determined 
conditions. 

Bis-intercalators like YOYO and ditercalinium, for in¬ 
stance, are molecules which have two intercalating por¬ 
tions linked by another chemical structure which some¬ 
times may interact with the DNA groove s 145 : 158 ! 159 . 
Actinomycin D is another example of a drug with 
distinct portions that interact to DNA by different 
modes, in this case including minor groove binding and 
intercalations^.^. 

Psoralen is an example of a drug which the binding 
mode depends on an external condition (sample illumina¬ 
tion). As explained before, the drug initially intercalates 
in DNA, but forms covalent bonds with the thymines if 
the sample is illuminated with UVA light. 

Hoechst 33258 is a fluorescent stain that can bind to 
DNA by intercalation or groove binding, with two dif¬ 
ferent sets of physicochemical parameters. In this case 
the drug concentration is the factor that determines the 
dominant binding modeSS. Some authors report a similar 
behavior for the intercalator doxorubicin, have founding 
a possibility of groove binding at AT-rich regions^. 


VI. CONNECTING MECHANICS TO PHYSICAL 
CHEMISTRY 

In this section we introduce the main subject of this re¬ 
view, the approaches developed to establish connections 
between the mechanical properties and physicochemical 
properties of DNA-ligand complexes. As stated before, 
the advantage in establishing this type of connection is 
the possibility to deduce one or more properties of a 
certain type (physicochemical properties, for example) 
knowing only the behavior of a property of the other 
type (the persistence or contour length, for example). 
With such connection(s), one can considerably reduce 
the number of different experimental techniques neces¬ 
sary to perform a robust characterization of the DNA 
interaction(s) with a certain type of ligand. This fact 
thus reduces the time and cost required for getting data, 
since less different equipments are needed and the num¬ 
ber of experiments that must be conducted can be con¬ 
siderably reduced. Furthermore, and perhaps more im¬ 
portant, the approach opens the possibility of comparing 
data obtained by means of very different experimental 
techniques, increasing confidence in the results. 

In single molecule stretching experiments performed by 
optical or magnetic tweezers, the typical result obtained 
is the force x extension curve of the molecule, from where 
the mechanical properties can be extracted by fitting an 
appropriate model (for DNA, the WLC model). We will 
show that if one knows how the contour and/or the per¬ 
sistence length varies as a function of the total concen¬ 
tration of ligand in solution (Ct) (which is the amount 
of ligand added in sample preparation), it is possible to 
deduce physicochemical properties such as the equilib¬ 
rium constants, the cooperativity degree, the exclusion 
number, etc. 


A. Historical survey 

Along the past years many groups have used 
single molecule techniques to identify the pos¬ 
sible binding mechanisms of DNA-ligand inter¬ 
actions and to extract physicochemical infor¬ 
mation of such interactions from these types of 
experiment s 56 ' 63 ' 65 ' 75 ' 90 -92,97 ' 98 ’ 137 ' 140 ’ 158,165 ~ 172 * 

To the best of our knowledge, the first attempt to 
connect mechanics to physical chemistry for DNA-ligand 
systems was performed in the early 80’s by the group 
of D. M. Crothers, who have measured the changes of 
the DNA contour length when interacting with various 
drugs (netropsin, distamycin, iremycin, daunomycin) as 
a function of the bound ligand fraction r, by using electric 
dichroism and a phase partition techniqu e 137 : 165 . Never¬ 
theless, they have not directly determined physicochem¬ 
ical properties from such data, a task which could only 
be performed with complementary analyzes and/or tech¬ 
niques. 

In 1996 Coury et has determined physicochemi- 
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cal properties from the contour length data of some DNA- 
ligand complexes, obtained using AFMi££. In fact, by 
determining the relative increase of the contour length of 
DNA complexes formed with intercalators such as dauno- 
mycin and ethidium bromide, the authors were capable 
to estimate binding parameters such as the equilibrium 
constant and the exclusion number. Similar approaches 
were used by Mihailovic et al. 767 and Rocha et al. 75 to 
extract physicochemical information of DNA-intercalator 
complexes by measuring the relative increase of the con¬ 
tour length, obtained using optical tweezers. 

Basically, the idea to perform such task is the following. 
Call A the natural distance between two DNA base pairs, 
which is ~ 0.34 nrn for B-DNA. When an intercalating 
molecule binds to this site, it increases such distance to a 
new value A + <5. Call Lq the bare DNA contour length 
and L the new length for a certain amount of bound 
ligand represented by the bound fraction r. One can 
promptly write the relation 


equilibrium dialysis, etc.) to evaluate such partitioning. 
There are, however, approaches that can be performed 
to bypass this problem, allowing one to use only single 
molecule stretching to characterize the interaction. In 
fact, one can estimate the bound ligand concentration 
from contour length changes if the length increase due to 
a single binding event (<5) is know n 166 : 167 . Alternatively, 
as a first-order approximation one can consider Cf ~ 
Ct in the binding isotherm if the DNA concentration 
in the sample is very low (because C b will also be very 
low in this case) 63-16 '. This approximation is much used 
in typical tweezers experiments that tether an individual 
DNA molecule and then rinse away any DNA molecules 
in solution prior to the introduction of a ligand. Such 
approach is convenient because it allows one to express 
the binding isotherm as a function of a directly accessible 
parameter (Ct), although it cannot be used always. A 
different approach was proposed originally by Rocha et 
al. in 2007—, which consists in manipulating Eqs. HM1E51 
and M to write the relation 


L = L 0 + N b S, (22) 

where N b is the number of bound ligand molecules. 

Observe that L 0 = N bp A, being N bp the number of 
DNA base-pairs. Therefore one can write the relative 
change of the contour length O as 


L - L 0 N b S 
L 0 N bp A 7r> 


(23) 


where 7 = 6/A. 

With Eq. [231 one can directly connect the mechanical 
parameter L to physicochemical parameters by express¬ 
ing the bound ligand fraction r by an adequate binding 
isotherm. In the case of intercalators, the more conve¬ 
nient binding isotherm is the Neighbor Exclusion Model 
(NEM) (see Section HID), since this isotherm captures in 
detail the neighbor exclusion effects which always follow 
intercalative binding. Nevertheless, there are two prob¬ 
lems that must be bypassed to use this approach. The 
first one is that in the NEM binding isotherm one can¬ 
not analytically isolate the parameter r to substitute in 
Eq. 1231 This problem, however, can be bypassed with 
numerical approaches, as will be discussed soon. The 
second and more serious problem is that, not only NEM, 
but all binding isotherms are written as functions of the 
free ligand concentration Cf, which is not a directly ac¬ 
cessible parameter. In fact, in general one knows only 
the total ligand concentration in solution Ct, the quan¬ 
tity used to prepare the sample, which is the sum of the 
bound and free ligand concentrations, i. e., 


C T = C f + C b . (24) 

The partitioning of Ct into C / and C b is not trivial to 
be measured and one usually needs other experimental 
techniques (microcalorimetry, absorption spectroscopy, 


C T = —B + 
7 


@(7 — n 0 + @) n_1 
Ki( 7 — nO) n 


(25) 


Such approach allows one to directly fit the contour 
length data without any approximation: one should just 
plot the total ligand concentration Ct in the y-axis and 
the relative increase of the contour length 0 in the a;-axis, 
such that Eq. [25] can used directly to fit the experimental 
data. In Fig. [2] we show an example of such fitting, 
performed originally in ref.' - for DNA complexes formed 
with the intercalator diaminobenzidine. Other examples 
can be found in refs . 75 and—2 for DNA complexes formed 
with the intercalators ethidium bromide and psoralen, 
respectively. 

One should note, however, that the contour length ap¬ 
proaches discussed above can only be used for intercala¬ 
tors. In fact, only intercalators increase the DNA con¬ 
tour length when bindin g 56 i 73 i 137 . An exception are lig¬ 
ands that facilitate or inhibit base pair formation, which 
can also be studied by contour length approaches simi¬ 
lar to those discussed above, using length changes as a 
marker of ligand binding^— The other common types 
of interactions between DNA and ligands, such as groove 
binding, electrostatic interaction or covalent binding in 
general do not affect the DNA contour length. In some 
cases, however, these kinds of interactions can cause DNA 
compaction with a decrease of the “apparent contour 
length” measured by force spectroscopy in the low-force 
regim e 65 ’ 97 : 98 . The concept of “apparent contour length” 
arises from the fact that, if the DNA molecule is partially 
compacted due to ligand binding, small forces in the en- 
tropic regime usually are not sufficient to fully stretch 
the molecule, and therefore the measured contour length 
will be smaller than the real one. Depending on the 
type of interaction, even high forces cannot be used to 
fully stretch the complexes and estimate the real contour 
length by fitting the WLC model^— The decrease of the 
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FIG. 2. Experimental result ( circles ) of Ct x 0 mesured for 
DNA complexes formed with the intercalator diaminobenzi- 
dine, and a fitting to Eq. [5S ^(dashed line). Observe that Eq. 
[25] fits well to the experimental data, returning the values of 
the physicochemical parameters N = 2.5 ± 0.6, Ki = (1.8 ± 
0.6) xlO 4 M -1 and 7 ~ 1. For this data Cb p = 2.4 fiM. 


“apparent contour length” upon increasing of the bound 
ligand concentration in general depends on intricate ef¬ 
fects such as the positional correlation of bound ligands. 
This fact makes it difficult to directly link the contour 
length data to a binding isotherm, although other kinds 
of analyses can be performed to study such interactions. 

On the other hand, the other basic mechanical 
property (the persistence length) is much more sen¬ 
sitive to other types of interactions, and in gen¬ 
eral changes for covalent bindin g 97 i 98 i 126 i 133 ; intercala¬ 
te bindin g 56 i 75 ~ 77 i 103 i 140 i 141 and groove/electrostatic 
bindin g 19 ' 56 ! 64 ' 65 ! 91 ' 92 ! 155 . This fact turns the persis¬ 
tence length into the ideal mechanical property to be 
choosen for monitoring DNA-ligand interactions and to 
be connected with the physical chemistry of such inter¬ 
actions. Nevertheless, such connection is not straight¬ 
forward as the one performed for the contour length of 
DNA-intercalator complexes. 

Only in 1998 the first attempt to connect the per¬ 
sistence length to physicochemical properties was per¬ 
formed by Rouzina and Bloomfield^, which can be syn¬ 
thesized in Eq. [2D] presented earlier. This model however 
was derived in the context of electrostatic interactions 
and attempt to explain the changes of the persistence 
length due to the negative charge neutralization in the 
DNA phosphate backbone*!— Recently, this model has 
been used to fit experimental data of DNA complexes 
formed with positively charged proteins such as HMG, 
HMGB1 and HMGB 9 91 ! 92 , with excellent agreement. 

The question now is: Can the changes of the persis¬ 
tence length be related to physicochemical parameters 
for any type of interaction? In the next section we dis¬ 
cuss an approach that can be used to perform such task. 


B. A general model to connect the persistence length to 
physical chemistry 

The DNA molecule partially covered by ligand 
molecules along its structure can be thought as an as¬ 
sociation of entropic springs in series. One type of 
spring is the bare DNA with its natural persistence 
length Aq, corresponding to the regions without bound 
ligands along the contour length of the molecule. The 
other type(s) of spring(s) is(are) the local complexe(s) 
formed between DNA and the bound ligand molecules. 
A simple phenomenological model to study the persis¬ 
tence length of DNA-ligand complexes that uses this 
assumption was proposed by Rocha—. Latter, it was 
rigourously demonstrated^ that a series association of 
n entropic springs with persistence lengths Aq, Ai, A 2 , 
..., A„_i results in an effective entropic spring with the 
effective persistence length Ae given by 


1 = fair) fijr) hjr) 
Ae Aq Ai A2 


(26) 


where fo(r), fi(r), f 2 {r), etc are specific functions of the 
bound ligand fraction r. 

The function fi(r) is in fact the probability of finding 
an entropic spring (a part of the DNA molecule) along 
the contour length with a local persistence length Aj&, 
which depends on the bound site fraction r 63 i 66 . 

In general, the approach proposed in Eq. [26] can be 
applied by following three steps: 

(a) One needs firstly to find the probability distribu¬ 
tion of the bound ligands, i. e., the set of functions fi{r). 

(b) The second step is to choose an adequate bind¬ 
ing isotherm that captures the physical chemistry of the 
system, and then plug such isotherm in Eq. [21] via the 
parameter r. 

(c) Finally, the third step is to use the equation con¬ 
structed in step (b) to fit the experimental data of the 
persistence length, extracting the physicochemical pa¬ 
rameters contained in the binding isotherm and the set 
of local persistence lengths Aj’s. 

To deduce the probability distribution mentioned in 
step (a), the easiest way is firstly identify how many dif¬ 
ferent entropic springs one needs in the model to correct 
reproduce the experimental behavior of the persistence 
length as a function of ligand concentration. The sim¬ 
plest behavior of this parameter reported in the litera¬ 
ture is a monotonic decay, found for example for the pro¬ 
teins HMG, HMGB1 and HMGB 2 91 ! 92 and for the drug 
cisplati n 97 ! 98 . This relatively simple behavior of the per¬ 
sistence length can be explained with a model consisted 
only by two entropic springs, one representing the bare 
DNA (local persistence length Ao) and the other repre¬ 
senting the local structure formed by the ligand molecule 
bound to the DNA (local persistence length A\). 

From now on let us consider a “site” the place ef¬ 
fectively occupied by a single ligand molecule (or by a 
single bound cluster of molecules, in the cases in which 
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the ligands bind to DNA forming clusters due to high 
positive cooperativity). One should note that even for 
single-ligand binding the sites are usually higher than 
one DNA base-pair, due to ligand size and/or neighbor- 
exclusion effects. A model with only two different types 
of entropic springs, as proposed in the last paragraph, is 
a one-site quenched disorder statistical model , since the 
probability distribution depends only on the occupancy 
of single sites along the double-helix, i. e., it does not 
depend on the correlation with the occupancy of nearest 
neighbor sites. 

Consider now a particular site choose randomly along 
the DNA. The probability of this site to be occupied by 
a ligand molecule is x = r/r m ax , with a local persistence 
length Ai] and the probability of this site to be unoccu¬ 
pied is 1 - x, with a local persistence length A o 63,66 . The 
effective persistence length can then be written as 


Ae A 0 A\ 


(27) 


and x = r/r ma x can be directly connected to a binding 
isotherm. 

Equation [57] was recently used by Crisafuli et al. to 
determine the physicochemical parameters of the DNA- 
cisplatin interaction from the persistence length data of 
these complexes2£i2£. One should observe that, since the 
exclusion number N is related to the saturated bound 
ligand fraction r max by N = l/r max , the electrostatic 
model proposed by Rouzina and Bloomfield^ (Eq. I5U1) 
is a particular case of Eq. [27] (it is just a matter of redefin¬ 
ing the physical interpretation of the constants Aj’s). 

There are other types of ligands that can induce 
a more intricate non-monotonic behavior for the per¬ 
sistence length as a function of the ligand concentra¬ 
tion. Probably the most known example is the bacte¬ 
rial protein HUiS, but some drugs such as catinonic 
cyclodextrins^, actinomycin D^i and hoechst 33258^ 
also induce such behavior. To account for the persistence 
length changes of the DNA complexes formed with these 
compounds, the one-site model discussed above does not 
work, and one needs to introduce at least one more en¬ 
tropic spring with other local persistence length (A 2 ), 
i. e., one needs a two-sites quenched disorder statisti¬ 
cal model , in which one must consider the probabilities 
associated with the occupancy of two nearest sites. In 
the context of a two sites model, there are therefore the 
following probabilities associated to the local persistence 
lengths: (a) two nearest sites unoccupied have local per¬ 
sistence length A 0 and probability Pq = (1-x) 2 . (b) Two 
nearest sites simultaneous occupied have local persistence 
length A 2 and probability P 2 = x 2 . (c) Finally, one site 
unoccupied and the neighbor occupied have local per¬ 
sistence length Ai and probability P\ = 1 - Pq - P 2 = 
2x(\-x). The effective persistence length can therefore 
be written as 


1 (1 — x) 2 2 x(l — x) x 2 

Ae Aq Ai A 2 ’ 


(28) 


and x = r/r max can be connected to a binding isotherm 
as usual. 

A last issue must be solved to complete the problem, 
both for monotonic and non-monotonic behaviors of the 
persistence length: one must write the binding isotherm 
as a function of a directly accessible parameter instead 
of Cf, as discussed in Section VIA, in order to eliminate 
the dependence in using other experimental techniques to 
estimate the ligand partitioning between the DNA (C&) 
and the solution (Cf). Although we have discussed some 
approaches to perform this task for intercalators in Sec¬ 
tion VIA, it is clear that a general approach is needed in 
order to contemplate the order types of ligands. 

In 2012 Siman et al. have firstly proposed a simple 
iterative solution of the binding isotherm^, which was 
promptly generalized by Cesconetto et al. in 2013 with 
the following method. Firstly choose a particular binding 
isotherm, for example, the Hill binding isotherm (Eq. 
[Til) . One can plug the relations x = r/r ma x and Cf = 
Ct - rCbp = Ct - rmax CbpX in this binding isotherm to 
write 


[Kj(C T - r max C bp x)] n 
1 + [. Ki(Ct - r m axCb P x)] n 

Observe that this equation can be solved numerically 
for known values of the constants, returning x for each 
value of Ct- Therefore, one needs to write a simple al¬ 
gorithm that uses a subroutine to solve Eq. [55] for initial 
guessed values of the constants, and uses the results re¬ 
turned for x plugged into Eq. [55] or Eq. [57] to fit the 
experimental data of the persistence length A as a func¬ 
tion of Ct, by using least squares fitting. With this ap¬ 
proach the problem is completely solved. Observe that 
any binding isotherm can be used to get an equation sim¬ 
ilar to Eq. [55] i. e., one needs only to choose a plausible 
binding isotherm that captures the physical chemistry of 
the interaction. 

Below we revisit some results recently obtained with 
this approach, showing that in principle it can be used 
to study any type of interaction. The only requisite is 
that such interaction changes the DNA persistence length 
as the ligand binds. All the experimental data were ob¬ 
tained by single molecule stretching performed with opti¬ 
cal tweezers in the entropic regime, with fittings similar 
to that shown in Fig. |T| (except those of Fig. [5] - see 
ref/ 1 '-). We also discuss the main features of the physics 
and chemistry of the interactions revisited here. A com¬ 
plete detailed discussion can be found in the original ar¬ 
ticles (and the references therein). It is worth to em¬ 
phasize that all optical tweezers measurements were per¬ 
formed in chemical equilibrium, waiting sufficient time 
for ligand equilibration before performing the stretching 






Extracting physical chemistry from mechanics: a new approach... 


13 


experiments. In addition, these measurements were per¬ 
formed with low forces (< 2 pN) and pulling rates (~ 0.1 
/im/s) in order to guarantee that the chemical equilib¬ 
rium is not affected by the stretching forces. 

In Fig. [3] we show the experimental data ( circles ) of 
the persistence length of DNA-cisplatin complexes as a 
function of drug total concentration in the sample Ct- 
Observe that, for convenience to fit with Eq. [271 we have 
plotted the inverse of the persistence length in this figure 
and in all subsequent ones. The fitting with the model 
(Eq. [27l) is also shown ( red solid line). In this case we 
have used the Hill binding isotherm (Eq. 1141) to perforin 
the fitting, extracting the physicochemical parameters A',; 
= (1.6 ± 0.2) x 10 4 M n = 3.6 ± 0.4, r max = 0.56 ± 
0.06 and A\ = (24 ± 4) nm. These results agree very well 
to those presented in ref.—, which were obtained using 
another fitting strategy, and as well as to results obtained 
from other experimental techniques 176,177 . In particular, 
the Hill exponent n = 3.6 indicates that cisplatin presents 
positive cooperativity in its interaction with DNA. 

Cisplatin and its analogues carboplatin and oxaliplatin 
form one of the most important class os compounds used 
in cancer chemotherapies, especially to treat head, neck, 
testicular, ovarian and non-small cell lung cancersiS. 
In aqueous solution, two chloride ions dissociate from 
the compound, followed by incorporation of two water 
molecules. This is the active state of the drug, which 
can bind to DNAiS. Many aspects of the DNA-cisplatin 
interaction are currently well established in the litera¬ 
ture, such as the mechanism of action of the compound 
as an anticancer drug, which consists in damaging the 
DNA molecule with adducts that form interstrand and 
intrastrand crosslinks^. These crosslinks hinder DNA 
replication by introducing strong structural perturba¬ 
tions on the double-helix such as bendings, partial un¬ 
winding and loop o 126 d 27 i 180 . These structural pertur¬ 
bations are closely related to the result found for the 
Hill exponent (n ~ 3.6). In fact, a positive coopera¬ 
tivity could be expected in DNA-cisplatin interaction, 
since the crosslinks and loops induced in the DNA by 
the drug approximate different strand segments as the 
drug concentration is increased, therefore increasing the 
probability of forming even more crosslinks and loops as 
cisplatin bind o 97 i 98 . A nearly similar mechanism was re¬ 
cently observed for the H-NS binding protein by Dame 
et al. , which have shown that a cooperative behavior in 
this case arises as an intrinsic property of DNA bridging 
due to duplex proximity!^— 

In Fig. [4] we show the experimental data ( circles ) for 
the inverse of the persistence length of DNA complexes 
formed with a monovalent cationic /3-cyclodextrin (6- 
monodeoxy-6-monoamine-/3-cyclodextrin) as a function 
of drug total concentration in the sample Ct, firstly pre¬ 
sented in ref—Cyclodextrins (CDs) are cyclic oligosac¬ 
charides composed of D-glucose units joined by glucosidic 
linkages. The /3 subtype consists of seven units and has a 
structure that resembles a truncated cone, with hydroxyl 
groups localized at the outer surface of the cone. That 



FIG. 3. Circles: inverse of the persistence length of DNA- 
cisplatin complexes measured by single molecule stretching 
experiments. Red solid line: a fitting to the model (Eq. 1271) 
using the Hill binding isotherm (Eq. 1141) . From this fitting 
we have found the physicochemical parameters K t = (1.6 ± 
0.2) x 10 4 M - 1 , n = 3.6 ± 0.4, r max = 0.56 ± 0.06 and A i 
= (24 ± 4) nm. For this data Cb p = 8.9 /iM. 



FIG. 4. Circles: inverse of the persistence length of DNA- 
cyclodextrin complexes measured by single molecule stretch¬ 
ing experiments. Red solid line: a fitting to the model (Eq. 
l28ll using the Hill binding isotherm (Eq. fl4l) . From this fit¬ 
ting we have found the physicochemical parameters Ki = (9 
± 1) x 10 4 M n = 3.7 ± 0.4, Ai = 8.4 ± 1 nm and A 2 = 
149 ± 21 nm. For this data Cbp = 11 /rM. 

gives CDs the property to be water-soluble and to have 
a relatively hydrophobic inner cavity able to partially or 
entirely accommodate polymers forming host-guest inclu¬ 
sion complexes^. Monovalent cationic /3-cyclodextrin 
is usually obtained by substituting one of the hydroxyl 
groups by an amino group. This molecule has been used 
for condensing DNA and introducing it into small vesicles 
for gene therapy applications^— 

Observe in Fig. [I] that for the complexes formed be¬ 
tween DNA and cationic /3-cyclodextrin the persistence 
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FIG. 5. Circles: inverse of the persistence length of DNA- 
HU complexes measured by single molecule stretching exper¬ 
iments (experimental data by van Noort et al —i, error bars 
are not available in this case). Red solid line: a fitting to 
the model (Eq. l28l) using the Hill binding isotherm (Eq. H4l) . 
From this fitting we have found the physicochemical param¬ 
eters Ki = (3.4 ± 0.4) x 10 7 M n = 3.6 ± 0.3, Hi = 
(8.5 ± 1) nm and Ai = (115 ± 13) nm. For this data Cb p is 
unknown due to the sample preparation procedure— It was 
left as an adjustable parameter, and the fitting returns Cbp ~ 
110 nM. 


length exhibits a non-monotonic behavior, and therefore 
we have used Eq. [2H to fit the data, together with the 
Hill binding isotherm ( red solid line). We have found the 
results Ki = (9 ± 1) x 10 4 M ~ 1 ,n = 3.7 ± 0.4, Hi = 
(8.4 ± 1) nm and A 2 = (149 ± 21) nm. The parameter 
fmax = 0.67 was known for this ligand such that we have 
fixed its value in the fitting— The value obtained for the 
Hill exponent n again indicates that the system is posi¬ 
tively cooperative, in this case forming bound clusters of 
~ 4 drug molecules at the binding sites— 

The results obtained for the bacterial protein HU in 
ref* 63 are somewhat similar, as shown in Fig. [5] We 
have used again Eq. [28] and the Hill binding isotherm to 
perform the fitting ( red solid line), and the experimental 
data ( circles) were obtained by van Noort et al. for this 
ligandi— From the fitting we have obtained the results 
K, = (3.4 ± 0.4) x 10 7 M - 1 , n = 3.6 ± 0.3, Hi = (8.5 
± 1) nm and A 2 = (115 ± 13) nm. Here again, r max 
= 0.11 is a known parameter and was maintained fixed 
in the fitting— For this data Cbp is unknown due to 
the sample preparation procedure— It was left as an 
adjustable parameter, and the fitting returns Cbp ~ 110 
nM. Here the fact that the persistence length increases for 
high protein concentrations agrees with results obtained 
in AFM images, which have shown the formation of rigid 
filaments^— 

At this point it is necessary to reflect on the use of 
the Hill binding isotherm in the analysis of DNA-ligand 
systems. The fact that we have found a Hill exponent 
n 1 for cisplatin, cyclodextrin and HU strongly in¬ 


dicates that relevant positive cooperativity is present in 
such systems. In fact, binding isotherms with no cooper¬ 
ativity such as the Scatchard model (Eq. IT21) or the basic 
Neighbor Exclusion Model (NEM) (Eq. fl5l) do not work 
in performing these fittings. The cooperative version of 
the neighbor exclusion model (Eq. 1161) in principle could 
be used, but we were not successful in performing the 
fitting in an easy way, founding numerical problems in 
solving the equation analogue to Eq. [29] for this bind¬ 
ing isotherm. In fact, the intricacy of Eq. 1161 somewhat 
limits its applicability in the fitting approaches discussed 
in this review. For this reason we use the Hill binding 
isotherm, a much simpler equation that also takes into 
account cooperativity effects. 

The first example of a non-cooperative system stud¬ 
ied with our approach are the DNA complexes formed 
with the drug Actinomycin D (ActD), firstly presented 
in ref—. This drug is a DNA ligand clinically used as 
an antibiotic and to treat some highly malignant can¬ 
cers, such as gestational trophoblastic disease^, Wilms’ 
tumord££ and rhabdomyosarcoma—i. The drug exhibits 
a complex interaction with double-strand DNA, present¬ 
ing two distinct parts which bind to DNA by different 
modes: while the phenoxazone ring intercalates, prefer¬ 
entially at the CG base pairs, the cyclic pentapeptide 
chains bind to the minor groove, usually forming hydro¬ 
gen bonds with the guanine bases^~-i— 

Here we clearly have the option of choosing different 
binding isotherms to perform the fitting. This fact il¬ 
lustrates the versatility of our approach, which returns 
consistent results even for different binding isotherms: it 
is required only to choose one that captures the basic 
physical chemistry of the system. In fact, if the DNA- 
ActD interaction is non-cooperative^, one can choose the 
Scatchard model or the basic (non-cooperative) neighbor 
exclusion model. Nevertheless, instead of the first option 
(Scatchard), we have chosen the Hill binding isotherm 
to perform the fitting. If everything is right, one should 
find a Hill exponent near unity (n ~ 1), since in this 
case the Hill model is just equivalent to the Scatchard 
one. Figure [ 6 ] shows the experimental data points ( cir¬ 
cles) and the fittings to Eq. [25] with the Hill model (red 
solid line) and with the neighbor exclusion model (blue 
solid line). From the first fitting (Hill), we find Ki = (1.5 
± 0.4) x 10 6 M - 1 , n = 1.1 ± 0.2, r max = 0.11 ± 0.01, 
Hi = (15.2 ± 0.6) nm and H 2 = (64 ± 25) nm. From 
the second fitting (NEM) we find Ki = (4.6 ± 0.5) x 
10 6 M " 1 , N = 4 ± 0.5 (the exclusion number for each 
bound ActD), Hi = (14 ± 2) nm and H 2 = (140 ± 16) 
nm. Observe that both fittings explain well the behavior 
of the experimental data. The results returned for the 
physicochemical parameters, although somewhat depen¬ 
dent on the chosen binding isotherm, are realist. The 
relatively high variability on the values found for some of 
these parameters is compatible to the variability found 
when using different experimental technique s 53 ’ 187 H ^ 2 . 

A relevant question that can be raised at this point 
is about the accuracy of our approach to treat systems 
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FIG. 6 . Circles: inverse of the persistence length of DNA- 
ActD complexes measured by single molecule stretching ex¬ 
periments. Red solid line: a fitting to the model (Eq. 1281) 
using the Hill binding isotherm (Eq. 1141) . From this fitting 
we have found the physicochemical parameters Ki = (1.5 ± 
0.4) x 10 6 M n as 1.1 ± 0.2, r max = 0.11 ± 0.01, Ai = 
(15.2 ± 0.6) nm and A 2 = (64 ± 25) nm. Blue solid line: a 
fitting to the model (Eq. 1281) using the NEM binding isotherm 
(Eq. 1151) . From this fitting we have found the physicochem¬ 
ical parameters Ki = (4.6 ± 0.5) x 10 6 M -1 , N = 4 ± 0.5 
(the exclusion number for each bound ActD), Ai = (14 ± 2) 
nm and A 2 = (140 ± 16) nm. For this data Cb P = 10.6 pM. 


with more than one binding mode, i. e., with two or 
more different sets of physicochemical parameters. Many 
compounds interact with DNA in this way, and recently 
we have successfully applied our fitting approach to the 
fluorescent dye Hoechst 33258^. 

The Hoechst stains, also known as bis-benzimides, are 
a family of fluorescent dyes largely employed to stain the 
DNA molecule in molecular biology applications, allow¬ 
ing one to visualize DNA with fluorescence microscopy. 
In addition, these compounds can be potentially used as 
anticancer drugs^ 2 £, since their strong interaction with 
DNA can impede the replication of the molecule. Many 
experimental techniques were employed over the past 
years to study the effects of the Hoechst 33258 sub- 
type on the DNA molecule. In particular, it was found 
that the ligand binds preferentially to the DNA mi¬ 
nor groove, especially at AT-rich region o 191 ' 192 .. Nev¬ 
ertheless, some authors have proposed that the ligand 
presents more than one binding mode to double-strand 
(ds) DNA-i 22 r-i 2 i, indicating the possibility of intercala¬ 
tion at GC-rich regions^ 92 ] 195 . 

With our fitting approach we were able to decouple 
the two main binding modes that Hoechst 33258 exhibits 
with DNA, by using a binding isotherm expressed as a 
sum of two Hill processes. We have determined the two 
complete sets of physicochemical parameters for each of 
the binding modes. In particular, we have found that 
the first binding mode (intercalation) is non-cooperative, 
with a Hill exponent ~ 1, while the second mode (groove 



FIG. 7. Circles: inverse of the persistence length of DNA- 
hoechst complexes measured by single molecule stretching ex¬ 
periments. Red solid line: a fitting to the model (Eq. 1281) 
using a sum of two Hill processes as the binding isotherm. 
From this fitting we decouple the two binding modes and find 
the physicochemical parameters K\ = (1.8 ± 0.4) x 10 6 M 
m = 1.1 ± 0.3, K 2 = (2.4 ± 0.2) x 10 5 M n 2 = 7 ± 
3, A\ = (300 ± 100) nm and r m ax = 0.32 ± 0.02. For this 
data Cbp = 20 /rM. 


binding) is highly positively cooperative, with a Hill ex¬ 
ponent ~ 7. Such conclusion is in agreement with previ¬ 
ous studies performed by other techniques (equilibrium 
dialysis and absorption spectroscopy) 193 i 194 . The two 
binding modes coexist in the entire concentration range 
studied here, but intercalation is dominant for Ct < 3 
/rM while groove binding is dominant for higher concen¬ 
trations. Figure [7] shows the experimental data ( circles) 
and the fitting ( red solid line). We have found that, for 
the intercalative binding mode, A'i = (1.8 ± 0.4) x 10 6 
M -1 , ni = 1.1 ± 0.3. On the other hand, for minor 
groove binding, we found I\ 2 = (2.4 ± 0.2) x 10 5 M 
-1 , n 2 = 7 ± 3. Also, we found A\ = (300 ± 100) nm 
and r ma x = 0.32 ± 0.02, which are global parameters 
independent of the binding mode. The parameter A 2 = 
28.3 nm was maintained fixed in the fitting since it is 
the saturation value of the persistence length, which can 
be directly determined from the data of Fig. [7] in this 
case. The error bars of the parameters obtained in this 
fitting are somewhat higher than those obtained for the 
other DNA-ligand systems presented above. This fact 
is due to the excess of adjustable parameters used in the 
fitting procedure in this case, because we have two differ¬ 
ent binding modes and consequently two sets of binding 
parameters. 

Finally we show an example of how our fitting ap¬ 
proach can also be used to analyze the contour length 
data of DNA complexes formed with ligands. GelRed is 
a fluorescent nucleic acid stain designed with the purpose 
of replacing the highly toxic ethidium bromide (EtBr) 
in gel electrophoresis and other experimental techniques 
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which depends on the fluorescence of stained DNA. When 
bound to DNA, GelRed has the same absorption and 
emission spectra of EtBr and, according to its manu¬ 
facturer (Biotium Inc., Hayward, CA, USA), it has the 
advantage of being much less toxic and mutageni c 196 ! 197 . 

Figure [S] shows the experimental data of the relative 
increase of the contour length 0 = (L - Lq)/Lq ( cir¬ 
cles ), firstly presented in ref. 7 —, and two fittings per¬ 
formed with Eq. [23] and two different binding isotherms: 
Scatchard ( red solid line) and NEM ( blue solid line). 
The two fittings are similar, returning equivalent physic¬ 
ochemical parameters and allowing one to conclude that 
GelRed interacts with DNA by bis-intercalation' 7 . For 
the Scatchard fitting, we have found K, = (1.8 ± 0.4) x 
10' M _1 , r ma x = 0.22 ± 0.03 and 7 = 2.2 ± 0.1. For the 
NEM fitting, we found AT t = (1.8 ± 0.3) x 10 7 M N 
= 3.7 ± 0.4 and 7 = 1.9 ± 0.1. The results obtained for 
this ligand with our fitting approach lead us to conclude 
that the GelRed dye is a bis-intercalator. In fact, the ex¬ 
clusion parameter N = 1 /r max indicates that each bound 
GelRed molecule effectively occupies ~ 4 DNA base- 
pairs, a value considerably higher than the results found 
for most monointercalators, and approximately twice the 
result for EtBr (which is ~ 2^ —). The equilibrium asso¬ 
ciation constant A',; is also higher than the result obtained 
for typical monointercalators (~ 10 5 M _1 ) 73 -' 5 - 147 ., an d 
within the range found for most bis-intercalators (10' to 
10 9 M" 1 ) 145 ! 198 - 201 . Finally, the result 7 ~ 2 is approx¬ 
imately twice the value obtained for typical monointer¬ 
calators, suggesting that each bound GelRed molecule 
increases the DNA contour length by ~ 0.68 nm, a re¬ 
sult also compatible to typical bis-intercalator s 198 ’ 200 . 
Observe that the bis-intercalators should increase ap¬ 
proximately twice the DNA contour length per bound 
molecule, since each ligand molecule contains two inter¬ 
calating portions. 

In summary, we have presented many examples of 
DNA-ligand systems analyzed with the proposed fitting 
approach. All the results obtained for the physicochemi¬ 
cal parameters are consistent with most studies found in 
the literature that have used many different experimen¬ 
tal techniques, from crystallography to fluorescence res¬ 
onance energy transfe r 53 i 176 i 187 ~ 189 i 193 i 194 i 202 ~ ^i. Thus, 
our fitting approach allows a direct comparison be¬ 
tween the results obtained from single molecule stretch¬ 
ing experiments to those obtained from typical ensemble¬ 
averaging techniques, which are usually used to char¬ 
acterize the physical chemistry of DNA-ligand interac¬ 
tions. A weakness of the presented approach that can be 
pointed concerns the error bars of the physicochemical 
parameters obtained from the fitting procedure, which 
can be a bit high if the number of adjustable parameters 
used in the fitting is high, as in the case of multiple bind¬ 
ing modes. Nevertheless in these cases one can maintain 
fixed some parameters previously measured with other 
techniques and perform the fitting using only the ad¬ 
justable parameters of interest. In this way, the fitting 
approach can be used to compare and verify results ob- 



FIG. 8. Circles: experimental data of the relative increase 
of the contour length 0 = (L - Lo)/Lo of DNA-GelRed com¬ 
plexes. Red solid line: a fitting to Eq. [23]using the Scatchard 
binding isotherm (Eq. 1131) . from which we have found the 
physicochemical parameters Ki = (1.8 ± 0.4) x 10' M _1 , 
r m ax = 0.22 ± 0.03 and 7 = 2.2 ± 0.1. Blue solid line: a 
fitting to Eq. 1231 using the NEM binding isotherm (Eq. 1151) . 
from which we have found the physicochemical parameters Ki 
= (1.8 ± 0.3) x 10 7 M~ 1 ,N = 3.7 ± 0.4 and 7 = 1.9 ± 0.1. 
For this data Cb p = 2.4 /rM. 

tained from very different experimental techniques, and 
can still be useful in the investigation of DNA-ligand in¬ 
teractions. 


VII. CONCLUSIONS 

We have reviewed important topics of the field “DNA- 
ligand interactions”, from DNA mechanics to DNA- 
ligand physical chemistry, emphasizing how one can con¬ 
nect the changes of the mechanical properties of DNA 
induced by the binding ligand to the physicochemical in¬ 
formation of such interaction. This type of connection 
is extremely relevant because it allows one to perform a 
robust characterization of the interaction both from the 
point of view of the mechanical properties and of the 
physical chemistry of the interaction by using only one 
experimental technique: single molecule stretching ex¬ 
periments. Moreover, the possibility of performing such 
connection reduces the time and cost required for getting 
results about a DNA-ligand system, since less different 
equipments are required and the number of experiments 
that must be conducted can be considerably reduced. 
Furthermore, and more important, it opens the possibil¬ 
ity of comparing the results obtained by means of very 
different experimental techniques, in special when com¬ 
paring single molecule techniques to ensemble-averaging 
techniques. 

In particular, we reviewed a fitting approach recently 
proposed by our group to connect the persistence length 
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of the DNA-ligand complexes to the physical chemistry 
of the interaction. Such approach in principle can be 
used for any type of ligand, from drugs to proteins, even 
if there are multiple binding modes. However, a test with 
sequence-specific ligands^ is still needed. In any case, 
the only requisite to try the approach is that the interac¬ 
tion must change the DNA persistence length as the lig¬ 
and binds, which usually occurs for all types of common 
interactions (intercalation, covalent binding, electrostatic 
driven interactions and groove binding) at least for some 
ligand concentration range. 
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